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Abstract. Following the general numerical analysis of Melita and Woolfson (1996), I showed 
in a recent paper that a restricted, planar, circular planetary system consisting of Sun, 
[T I ■ Jupiter and Saturn would be captured in a near (2:1) resonance when one would allow 
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for frictional dissipation due to interplanetary medium (Haghighipour, 1998). In order to 



>■ ■ analytically explain this resonance phenomenon, the method of partial averaging near a 

^ ■ resonance was utilized and the dynamics of the first-order partially averaged system at 

Q ■ resonance was studied. Although in this manner, the finding that resonance lock occurs 

0\ • for all initial relative positions of Jupiter and Saturn was confirmed, the first-order par- 



p i' tially averaged system at resonance did not provide a complete picture of the evolutionary 

dynamics of the system and the similarity between the dynamical behavior of the averaged 



O 
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d ■ system and the main planetary system held only for short time intervals. To overcome 



these limitations, the method of partial averaging near a resonance is extended to the 
second order of perturbation in this paper and a complete picture of dynamical behavior 
of the system at resonance is presented. I show in this study that the dynamics of the 
second-order partially averaged system at resonance resembles the dynamical evolution of 
the main system during the resonance lock in general, and I present analytical explanations 
for the evolution of the orbital elements of the main system while captured in resonance. 

Key words: planetary dynamics , resonance capture , averaging . 
*E-mail:nader@lula.physics. missouri.edu 



1. Introduction 

Inspired by ttie results of extensive numerical investigations by Melita and Woolfson 
(1996) and in the framework of a restricted planar circular three-body system, I have 
recently shown that a planetary system consisting of a star and two planets, subject to 
dynamical friction due to a freely rotating uniform interplanetary medium, will be captured 
in resonance when the inner planet is more massive (Haghighipour, 1998, hereafter paper 
I). As an example of such a planetary system, the three-body system of Sun- Jupiter-Saturn 
was numerically integrated and it was shown that the orbital period of Saturn would be 
captured in a near (2:1) commensurability with that of Jupiter (Figure 1). In this system, 
the motion of Sun was neglected and Jupiter was assumed to be moving on a circular orbit 
with a known constant frequency. As a result of this resonance trapping, the eccentricity 
and semimajor axis of Saturn's osculating ellipse and also its orbital angular momentum 
and total energy became essentially constant (paper I). 
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Figure 1. Graph of the ratio of orbital period of Saturn to that of Jupiter versus time. The system is captured in 
a near (2:1) resonance. Integrations were performed on a timescale of 10 [T j /Iti) years where T j is the orbital 
period of Jupiter. The system was integrated with Sun at the origin of an inertial plane-polar coordinate system, 
Jupiter initially on the x-axis and Saturn initially at (a, e, 9, v) = (1.8380462, 0.0556, 0, 0) where a and e 
are the semimajor and the eccentricity of Saturn's osculating ellipse, respectively, and V is its true anomaly. The 
density of the interplanetary medium was taken to be equal to 16 times the mass of Jupiter uniformly spread in 
a spherical volume with 50 au radius and the masses of Sun, Jupiter and Saturn were taken to be constant and 
equal to their present values. 



This paper is devoted to the analysis of the dynamical behavior of the planetary system 
described above while captured in resonance. A preliminary analysis of the dynamics of 
this system was presented in paper I using the principle of averaging ; in fact in paper 
I, the first-order partially averaged system at resonance was obtained and its dynamics 
was studied. It turned out that in the first order of approximation, the partially averaged 
system was Hamiltonian with a periodic potential function that guaranteed the occurrence 
of resonance lock for all initial relative positions of the two planets. The dynamics of the 
first-order averaged system agreed with that of the main three-body system only for short 
time intervals and it was mentioned that in order to obtain an agreement which would 
hold for longer times, one had to extend the analysis to the study of the dynamics of the 
second-order partially averaged system at resonance. 

In this paper, a complete picture of the dynamical behavior of the three-body system 
of Sun-Jupiter-Saturn during resonance capture is presented. This is accomplished by ex- 
tending the calculations to the second-order partially averaged system at resonance, where 
one can also explain the behavior of Saturn's orbital elements and angular momentum 
during the resonance lock. 

The model under investigation is presented in section 2. Section 3 contains a brief 
discussion of the method of partial averaging near a resonance and in section 4, the appli- 
cation of this method to the Sun- Jupiter-Saturn system is discussed. Section 5 concludes 
this study by reviewing and summarizing the results. 

2. The Model 

The model under investigation is a restricted planar circular three-body system con- 
sisting of Sun, Jupiter and Saturn. In this model, the motion of Sun due to gravitational 
attraction of Jupiter and Saturn is neglected and Jupiter is assumed to have a uniform 
circular motion with a known constant period Tj. The only source of dissipation in this 
system is the dynamical friction due to interplanetary medium. It has been shown by 



Dodd and McCrea (1952) and also by Binney and Tremaine (1987) that the effect of the 
dynamical friction on the orbital motion of a general spherical body b with mass m^ around 
a general star with mass M appears as a deceleration with a magnitude given by 



R = ^2 ° ^b In 
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i + {^r , (1) 



rel 

and a direction that is opposite to the velocity of body h with respect to the medium. 
In this equation, Py is the uniform density of the interplanetary medium, Q is Newton's 
constant and S = rj^(m^/2M)^/"^, where r^, is the radial distance of b to the central star. 
The relative velocity Vrei in equation (1) is given by Vrei = V — Vm, where V is the velocity 
of the body and Vm is that of the interplanetary medium at the position of the body. For 
the planetary system presented in this paper, it is assumed that the interplanetary medium 
is freely rotating around the central star in the same sense as body h. Therefore, Vm is 
perpendicular to the position vector of h and has a magnitude equal to {GM/r^^Y''^. 

In an inertial coordinate system with its origin on Sun, the equation of motion of 
Saturn in dimensionless form is given by(paper I) 

iL^^L ,!fzlj) n . (2) 

at^ r'^ \r — r j\'^ 

In this equation, e = m j/M, where m is the mass of Jupiter, f j is the unit vector along 
Jupiter's orbital radius and TZ = TZVrei/^rei where 

^ = ^ ln(l + Br'Vii) , (3) 

is the dimensionless magnitude of Saturn's deceleration due to dynamical friction and Vrei 
is dimensionless magnitude of its velocity with respect to interplanetary medium and is 
given by 

Vrei = f^ + r\e - Umf, (4) 

where Um = r"^'^, is the dimensionless angular frequency of the medium at distance r 
from Sun. Denoting the mass of Saturn by rM^, one finds that A = 2Tv{p^^r^/M){mg/M) 
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and B = 2~^/^(mg/M)~'^/^. For any physical system A « 1 and d = AB « 1. It is 
important to mention that in deriving equation (2) a set of units has been chosen such that 
{Tj/27r)^ = r^ /QM and all positions and time variables have been scaled by the orbital 
radius of Jupiter (r^) and its orbital period divided by 27r, respectively (see paper I for 
details). 

3. Partially Averaged System at Resonance 

In order to apply the method of partial averaging near a resonance to the problem pre- 
sented in this paper, it is useful to write equation (2) in terms of the Delaunay action-angle 
variables (Appendix A). Denoting the semimajor axis and the eccentricity of Saturn's os- 
culating ellipse by a and e respectively, the appropriate Delaunay variables for the relative 
motion (2) are given by L = a^'^, G = L(l — e^)^' ^, I = u — e sinu and g = 9 — v, where 
u,v and / are the eccentric, the true and the mean anomalies of Saturn's orbit, respec- 
tively (Kovalevsky, 1967; Sternberg, 1969; Hagihara, 1972). In these relations, L and G are 
action variables and / and g are their corresponding angular variables. From Appendix 
A, it is evident that For the pure Kepler problem, i.e. in the absence of perturbations in 
equation (2), the solutions of the Delaunay equations are given by constant L,G and g 
and 

/ = uj^t + l^ , (5) 

where w^ = L~^ is the Keplerian frequency and /„ is a constant of integration. 

The problem under investigation here is a Kepler system with period T = l-n juOj^ 
that is perturbed by forces that involve a period 7^. At resonance, it is necessary, but 
not sufficient that T and 7^ become commensurate, i.e. relatively prime integers n and 
n' exist such that n7^ = n'T . To explain the dynamical behavior of the system near the 
resonance manifold in phase space, it is more convenient to imagine a system with a single 
degree of freedom and action-angle variables (X, G) such that 
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J- eP(J,e,t,e), (6) 

and 

e = 0(J) + eAr(J,e,t,e), (7) 

where V and M are explicitly time-dependent with period Tp. 

In the absence of perturbation (e = 0), equations (6) and (7) represent a Hamiltonian 
system whose action variable 1 has a constant value X^ . At resonance, this action variable 
fluctuates around its unperturbed value. Let us denote these fluctuations by 

J-J„=e'P, (8) 

where e' is a small parameter which has yet to be determined. In a similar fashion, let us 
show deviations of O from its unperturbed value Vt{X^^)t by 

e - o(jjt = $. (9) 

Substituting for X and Q from equations (8) and (9) in equations (6) and (7), the dynamical 
equations of the system will be given by 



e dV 

V = -V{I„e,t) +eV—{I^^,e,t) +0{ee'), (10) 

and 

/90 1 f)^0 

^ = '"^h^^o) + eAr(J„,e,t) + ^e''^'|^(^o) + 0{ee') + 0{e'') . (11) 

Using the averaging theorem (Sanders, 1985; Wiggins, 1996; Chicone et al., 1996a&b 
and 1997a&b) , we can now study the averaged dynamics of the system (10) and (11) at 
resonance by introducing the averaging integration 



1 r'^'^p 
Q = ^r Q{T,MWt + ^.t)dt. (12) 



for any function Q(X, G,t). Also, the averaging theorem requires V and $ to have the 
same power of perturbation parameter in the lowest order, that is e' = e . The partially 
averaged dynamics of the system (6) and (7) at resonance is now given by 

P = ei/2p(J,„$,t) +e^( — (J„,$,t)), (13) 



and 



dn 



1 ^9 / d^n 



^ = e'/^'^^{I,)+e[mi„^,t) + -V^(^^I„))\ , (14) 

where () represent an averaged quantity. Solutions to equations (13) and (14) are expected 
to represent the average dynamical behavior of the perturbation system (6) and (7) over 
time intervals of order e~^/^ (Sanders et al., 1985; Wiggins, 1996; Chicone et al., 1996a&b 
and 1997a&b). 

4. Second Order Partially Averaged System at Resonance 

Mathematically speaking, equation (2) represents a dynamical system with pertur- 
bation parameters e and d (see Appendix A). In a parameter space, e and d represent 
a plane where the full two-parameter perturbation problem (2) must be studied. This, 
however, requires consideration of all curves in this parameter space. However, in order to 
simplify the calculations and to avoid complexities, I consider a linear relation between e 
and 6. That is, 6 = eA, with A fixed. This will simplify the analysis by reducing the two- 
parameter perturbation problem (2) to a perturbation problem with only one parameter, 
e. 

At resonance, the orbital periods of Saturn and Jupiter become commensurate and 
therefore, the semimajor axis of Saturn's osculating ellipse becomes essentially constant. 
Following the formalism presented in section 3, L and / can be written as 

L = L^ + e^/^D , (15) 

and 

I = —t + ^ , (16) 





where D represents the deviation of L from its resonant value L^^ and ip represents the 
deviation of / from its resonant value L~^ t. Dynamics of the system near resonance is 
obtained by writing equation (2) in terms of the new action-angle variables {D, G, </?, (7, )• 
The result is (paper I) 

D ^ -e^/^Fii - eDFi2 + 0{e'^/'') , (17) 

^ = -£^22+ 0(£3/2) , (18) 

^ = -e^- (^) + - (^ + ^^32) + 0{e^'-) , (19) 

g =5^42+ 0(£=^/2) , (20) 

where 

F^^{L,,G,L-H + ^,g) = ^l^(L,,G,L;H + ip,g) - ATZUL^.G, L;H + cp, g) , (21) 

F22{L,,G,L;h + p^,g) = -^{L,,G,L-H + ^,g) - AnG{L,,G,L-H + ^,g) , (22) 

F^2{L,,G,L-H + p,,g) = -^(L,,G,L;H + ip,g) + Ani{L,,G, L;H + if, g) , (23) 
F^2{L,,G,L;'t + ^,g) = -^{L^,G,L;H + ^,g) + AngiL^,G,L;H + ^,g) , (24) 

and 

Fi2{L,,G,L;H + ^,g) = -^{L,,G,L;H + ^,g) . (25) 

In equations (21)-(25), the effect of gravitational attraction of Jupiter appears in the 
external Hamiltonian Tiext, where 

Hext = - rr^-n ■ (26) 

Also the contribution of the frictional force of the interplanetary medium is given by 
(Appendix C) 

Kl = |(l-e")-"2(K,yl-K„B) , (27) 



7Jo=i 



T^a; sin(gf + v) — IZy cos(q + v) 



(28) 



ni = -^a-^/^n^C + '}ZyW) , (29) 

^3 = ^ (^-^ + 7^y^) , (30) 

where 

A = sin {g + v) + e sin g , (31) 

B = cos {g + v) + e cos g , (32) 

C = (— 2e + cos t) + e cos^-O) cos {g + v) + {2 + e cos -0) sin (g + v) sin {i , (33) 

W = (— 2e + cos V + e cos^-O) sin {g + v) — (2 + e cos -0) cos {g + v) sin v , (34) 

/ sinO \ / N 

t = cos g + [ -— sm {g + v) , (35) 

V 1 + e cos V / 

J^ = sin5r-(— _ j cos(5r + 'f)) , (36) 

V 1 + e cos V / 

and T^x and 'Ky are the Cartesian components of 7?.. For the three-body system of Sun- 
Jupiter- Saturn in a (2:1) resonance, the dynamics of the second-order partiaUy averaged 
system at resonance is given by 

b = -e^/^Fii - eDFi2 , (37) 

G = -eF22 , (38) 

,^^^,m(^4)^,(^4^^P:^ , (39) 
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g = eF42 . (40) 



where 



r.47r 



F,^{G,g,^) = ^ [ F,j{G, l^t + ^,g,t)dt . (41) 

In order to study these equations, one needs to calculate Fij. From equations (21)- 
(25), calculations of these terms involve the effect of gravitational attraction of Jupiter 
and also the frictional effect of the interplanetary medium. Appendices B and C contain 
these calculations, respectively. From Appendix C and to the lowest order in eccentricity, 
the quantities IZl, T^g, T^i and IZg are given by 
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Ul^ --a'e' cos l{l-- cos^ /) ^'^ , (42) 

Ug - - - a e^ cos / (l - - cos^ /) ^^^ , (43) 

7^^ - - e^ (137 sin / - 23 sin 3/) (l - - cos^ /) ^^^ , (44) 

7^g = -Gae'^ (61 sin / + 23 sin 3/) (l - - cos^ /) ^^^ . (45) 
8 4 



Substituting for / from equation (16) and using the averaging integration (41), the averaged 
values of these quantities wiU be equal to zero. The non-vanishing contribution of friction 
appears as terms proportional to e^. Since at resonance the numerical value of e is only 
a few percent (paper I) , the dynamical friction will not play a vital role in the dynamical 
evolution of the averaged system. Therefore, in the rest of these calculations, I neglect the 
contribution of dynamical friction at resonance and will consider only the gravitational 
attraction of Jupiter as the main source of perturbation. Also, in order to simplify the 
calculations, all quantities will be expanded only to the first order in eccentricity. 

Neglecting the frictional perturbation of the medium, the quantities Fij will then be 
obtained only from the external Hamiltonian Tiext- From Appendix B these quantities are 
given by 



i^ii = 4^11 sin(2y: + g) , (46) 

% 

^2 = —^ ^11 sin (2v? + g) , (47) 



5/2 

etto 



-^22 



2a2 



4r + 3 4r 



sin{2ip + g) , (48) 



^32 = ^ ^11 cos (2v? + g) , (49) 

2eao 

G 

Fa2 = ^ — 3 CTii COS (2(/7 + g) , (50) 
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where 

CXD 

CTll = J2 

h=0 



,(3/2) 
'22 



a 



(5/2) 
22 



CXD 



^1=0 



ni + h) 


.<wr(f) 


"r(f + /i) 


[<wr(f) 


nl + h)' 



1 + 



2/1 + 3 



/i + l 



2/1 + 5 



4a2 V /i + 2 



1 + 



2h + 5\ /2/1 + 3 



1 



4a2 V /i + 2 
1 /2/1 + 5 



3 + 



/i + l 
2/1 + 7 



h + 2 



1 



3 2/1 + 9 



4a2 h + 3 



(51) 
(52) 
(53) 



a^^/i!r(|) r 4a2V/i + l 
Dynamical equations of the second-order averaged system at resonance can now be written 
as 



D = - 
G = - 



ee 



1/2 



^11 1 + 



e 



1/2 



2 1/2 



V' 



^ 



ee 

3 £1/2 



a, 



(3/2) 



22 



D + 



+ 3(J22 
6£ 



(5/2)- 



5/2 
tto 



D" 



£>) sin(2(/? + g) 

sin (2(/7 + g) , 
1 



12 e 



(Til cos(2(/? + (7) 



Ge 



(Til COS (2(/? + (7) 



(54) 
(55) 

(56) 
(57) 



2ea3 
where the overbars have been dropped for the sake of simphcity. 

Equations (54)- (57) allow us to analyze the averaged dynamics of the perturbation 
system (2) at resonance. Numerical integrations of these equations indicate that the rate 
of change of orbital angular momentum of Saturn (G) is almost zero (Figure 2) and its 
resonant value agrees with the value obtained from the direct integrations of the main 
system (paper I). ^ 0.2 0.4 0.6 




0.2 0.4 0.6 



Figure 2. Graph of orbital angular momentum of Saturn versus time at resonance. The initial conditions are 
given by {D, G, (/?, g) = (0, 1.275, -5.4, 0.882) and the timescale is 10'^{Tj/2tv) years. 
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One can also use G = L(l — e^)^'^ along with equation (15) to show that in the first 
order in eccentricity and second order in perturbation parameter e^'"^, 



de 



Ge 

2^ 



a 



(3/2) 



22 



+ 3a 



(5/2) 



22 



2Ga-^'''aii 



sin(2(/7 + g) ■ 



(58) 



Figure 3 shows that at resonance, the right hand side of equation (58) is essentially zero 
which implies that eccentricity becomes almost constant at resonance. 
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Figure 3. Graph of the rate of change of eccentricity of Saturn's osculating ellipse versus time. This graph shows 
that at resonance e is essentially constant. The initial conditions and the timescale are the same as Figure 2. 




Let us now turn our attention to evolution of action variable D with time. As men- 
tioned in section 2, D represents the deviation of L from its Keplerian value. One, therefore, 
expects L and D to have the same frequency. In a comparison with the result of the first 
order partially averaged system at resonance presented in paper I where the similarity 
between D and L was only pronounced for time intervals of order £~^' ^, the second-order 
partially averaged equations reveal this similarity in general (Figure 4). 
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Figure 4- Graph oi D (left) and the action variable L (right) versus time while the system is at resonance. The 
original system shows a periodic behavior (right) with a frequency which is in a very good agreement with the 
second-order averaged system at resonance (left) . The initial conditions for the graph of D (left) are the same 
as Figure 2 and for the graph of L are the same as Figure 1. The timescale for both graphs is 10 [Tj/2n) 
years. 

It is necessary to mention that the graph of L versus time shows two different fre- 
quencies (Figure 4). One is the frequency of its envelope and the other is the Keplerian 
frequency of Saturn. From the principle of averaging presented in section 3, it is the fre- 
quency of the envelope that has to be similar to that of D. Figure 5 shows the frequency 
of the envelope and also the Keplerian frequency of action variable L. The (2:1) commen- 
surability is simply obtained by calculating Saturn's Keplerian period from this figure and 
multiplying it by the timescale 10^(T^/27r). 
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Figure 5. Graph of the action variable L versus time while the system is at resonance. There are two frequencies 
associated with this quantity : A slow frequency related to the time variation of its envelope and a fast frequency 
which is equal to the Keplerian frequency of Saturn at resonance. The near (2:1) commensurability is simply 
obtained by dividing the time interval of this graph by the immber of revolutions of Saturn during this interval 
and multiply the result by 10 [Tj /Iti). 
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5. Summary 

It has been shown by Haghighipour (1998) and in a more extensive way by Mehta 
and Woolfson (1996) that the three-body system of Sun- Jupiter-Saturn wiU be captured 
in a near (2:1) resonance when one aUows for the frictional force of the interplanetary 
medium. In an attempt to study these results analytically, the method of partial averaging 
near a resonance was employed and the system was considered to be restricted, planar and 
circular. The first-order partially averaged system confirmed the occurrence of resonance 
capture for all initial relative positions of the two planets (Haghighipour, 1998); however, 
it was unable to fully illustrate the dynamical behavior of the system. The agreement 
between evolutionary dynamics of the first-order partially averaged system at resonance 
and the main system was observed only for time intervals of order {m ^ /M)~^''^. In order 
to obtain a better agreement, dynamics of the second-order partially averaged system at 
resonance has been studied. Numerical integrations of the averaged equations at this order 
have confirmed that the orbital angular momentum of Saturn and the eccentricity of its 
osculating ellipse become essentially constant at resonance and demonstrated the close 
similarity between the dynamics of the averaged system at resonance and the main system 
over long time intervals. 

The effect of dynamical friction has not appeared in this study. It has turned out 
that the first non-vanishing terms of the averaged frictional force will be proportional 
to e"^, which has been neglected in these calculations. It is important to mention that 
it is only through dynamical friction that one can study the effect of the density of the 
interplanetary medium on the dynamics of the system. However, calculations presented in 
this paper indicate that considering dynamical friction as the only source of perturbation, 
the density of the medium will not play a vital role in time variation of the orbital elements 
of Saturn at resonance. It is expected that accretion will provide the appropriate ground 
for the study of the effects of the density of the medium on resonance phenomena. 
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Appendix A 

Let us assume that the motions of Jupiter and Saturn are co-planar. In a plane-polar 
coordinate system with its origin on Sun, the equation of motion of Saturn, i.e. equation 
(2) can be written as (paper I) 



Pr = r 



Pr 

and 
Pe 

where 



r'e 



\r — r , 



[r - cosie-Oj)] - {TZ^cosO + 7^ySin^) 



— £ 



\r — r , 



s\n{e-ej) - r{-n^ sine + TZyCosO) 



7^, = -^ Ml + Br^Vt^i) [rcose - r{e - uj^)sine] , 

^rel 

TZy = -^ ln{l + Br^V^^i) [^sin6' + r(^ - c^^) cos6'] . 



V 



(Al) 
(A2) 
(A3) 

(A4) 

(A5) 
(A6) 



The Delaunay variables presented in section 4 are related to the plane-polar coordi- 
nates by Pq=G^ = r(l + ecos-O) ,Pr = e sin-y/G and g = 9 — v. From these equations, 
the dynamical equations (Al)-(A4) can be written as 



^ =a(l-e2)-i/2 
at 



FrC sin-f) -f F^) (1 + e cos-O) 



and 



dG 
dt 


= rFe , 






dl 
'dt ~ 


...+Ia-'/^ 


Fr 


dg 

dt 


1 

e 


a(l-e2)" 


1 


/2 



Fr(— 2e + cos-O + ecos^-O) — Fg (2 + e cos-O) sin-O 



-F^ cos-O + Fefl + -) 

\ 1 + e cosv/ 



sint) 



(A7) 

(A8) 
(A9) 

(AlO) 
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where 



Fr = e 



cos(^ - Oj 






sm{e-9j] 



\r — r , 



{IZx cos 6 + IZy sin 6) 



+ {IZx sin 6 — IZy cos 6) 



(All) 
(A12) 



and the Keplerian frequency of the osculating eUipse is given by o;^ = L~^ . The quantities 
Fr and Fg are associated with perturbations due to gravitational attraction of Jupiter and 
also dynamical friction of the interplanetary medium. 



Appendix B 



The external Hamiltonian Tiext can be written as 



Hext = -[r^ - 2rcos{9 - ^J + l] 



-1/2 



(Bl) 



Following the simplifications presented in paper I, quantities r and cos(6' — 9j) are approx- 
imately given by 

r ~ a ( 1 — e cos /) , (B2) 



and 



cos(6' — Oj) ~ cos{l + g — 6j) + e cos{2l + g — 6 j) — cos{g — 6 



where the O(e^) terms have been neglected. In this approximation, 



Hext^ -[l + a^-2acos{l + g-ej)] 



-1/2 



1 2acos/ + cos(2Z + (7 — 6*^) — 3cos((7 — 6*^ 

2 l + a2 -2acos(/ + 5f-6'^) 



(B3) 



(B4) 



From equations (21)-(25), contribution of gravitational attraction of Jupiter in dynamics 
of the second-order partially averaged system at resonance appears as partial derivatives 
of Hext with respect to L, /, G and g. Substituting a^ by L and e by (1 — G^/L^)^/^, one 
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will obtain 

dHext 



■a 



1/2 



1 + a^ -2acos(/ + fif-6',' 



-3/2 



Za cos / + — coslg — , 

e 2e 

1 



- 2 cos(/ + 5f - 6' J cos(2/ + g-Oj 



2e 



+ - ea 



3/2 



-5/2 



1 + a^ -2acos(/ + 5f-6'J 2(1 + 2a^) cos l-8acos{g-6 



- cos(3/ + 2£? - 29 j) + 3cos(/ + 2^? - 20,; 



(B5) 



Oriext 

di 


~a 


1 + a^ -2acos(/ + 5f-6',) 


-6/-2 


ea sin I -\- sirv{l -\- g — 9 j) 




+ esin(2/+^-0^) 




3 o r - T-5/2 




+ 2'" 


1 + a^ -2acos(/ + 5f-6'J 


sm{l + g-e,) 






2a cos / + cos(2/ + 5f - 6*^ ) - 3 cos(5f - 0j ) 



(B6) 



dn 



ext 



G 



dG 2e 



1 + a^ -2acos(/ + 5f-6',; 



-3/2 



2 a cos/ + cos(2/ + (/ — 6*^) — 3cos(5r — j] 



and 



dg 



'.a 



1 + 0? — 2a cos(/ + g — 9 j] 



-3/2 



- ea 
4 



(B7) 



(B8) 



sin(/ + 5f-6'j) + -e sin (2/ + ^f - 6'^ ) - 3sin(5f-6'j 

"I -5/2 

1 + a^ -2acos(/-f5f-6'j) 
4 sin / - 2a sm{g - 9 j) + 3 sin(/ + 2^? - 20 J 

- 2a sin(2/ + g-9j) - sin(3/ + 2g- 29 j 



Quantities above can be simplified noting that at resonance a^ ~ 1.625 (paper I). 
Therefore one can use the relation 



(1 -2^cosa + C2^ ^ = J2C^(^osa)C' 

q=0 



m < 1 



(B9) 
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with ^ = a ^ and a = {I + g — 9j), to write equations (B5)-(B8) in a simple form. In this 
relation, C^ (cos a) are Gegenbauer polynomials which are given by 

r^Xr X s^T{X + h)T{X + q-h) . -. 

C^icosa) = }_^ h\{q-h)\[V{\)]^ cos[{q-2h)a\ . (BIO) 

It is important to note that integration (41) applies certain restrictions to the admissi- 
ble values oiq and h. Due to the harmonic nature of the terms produced by equation (BIO), 
the only terms that survive the integration (41) are the ones proportional to cos{2l+g — 6 ,) 
and sin(2/ -^ g — 9 j). On the other hand, expansion of quantities (B5)-(B8) using Gegen- 
bauer polynomials, will reslut in producing terms proportional to cos[(q' — 2h){l + g — j)] 
and sm{{q — 2h){l + g — 6j)]. After simplifying equations (B5) to (B8), one realizes that the 
only terms that can have non-zero contributions are the ones with angular arguments equal 
to / , {g — 6 j) , {l + 2g — 29 j ) , (21 +g — 9 j) and {31 + 2g — 29 j ) which require q to be equal 
to (2/i±l) , {2h±2) , (2/i±3) , (2/i) and (2/i±l), respectively. Denoting the contributing 
part of the external Hamiltonian by 7Y*^^, the non- vanishing terms in equations (B5)-(B8) 
can be written as 

^ ^ - -^ ^11 cos(2/ + g-9,) , (Bll) 

'-'^ 2eao 

^^ ^ ±a^^sin{2l + g-9,) , (B12) 







and 



^ c. ^ [ag/^' + 3 aS/'>] M2l + 'J-9.) , (B13) 



Sni^t _, G £,jj cos(2i + 9-9,) , (B14) 



r-^ 



dG 2 e a3 







where an , cr22 and ^22 are given by equations (51)-(53). Substituting equations 
(B11)-(B14) in equations (21)-(25) and integrating the results using integration (41), 
-^117 -^127 -^22, -^32 and F42 can be written as in equation (46)-(50). 
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Appendix C 

I have shown in paper I that equations of motion in terms of the Delaunay variables 
can also be written as 



dL 

'dt 
dG 
'dt 

(U 
~dt 



- e —^ + sAUl 



dl 
dg 



+ eAUa 



1 d'hiext A ,^ 



and 



dg dHext , . ^ 



(CI) 
(C2) 

(C3) 
(C4) 



Comparing equations (C1)-(C4) with equations (A7)-(A10), 7^^,7^(3,7^^ and TZg can be 
written as in equations (27)-(30). After substituting for TZx and TZy from equations (A5) 
and (A6) in these equations, the effect of the frictional force of the interplanetary medium 
on dynamics of the averaged system appears as 



7^L 

ni-- 
n 



BVt 



^ (1 - e^)-^/^ [a (1 - e^) [6 - Um) + ersino] ln(l + Sr^V, 



^rel) 5 



r 



-^^(e-u:^)Hl + Br^Vtei) , 



(C5) 
(C6) 



ra 



rel 
-1/2 



eBV^i ^ 



r (— 2e+cos'0+ecos v)—r {6—uJm) C^^+ecosv) sin-O ln{l+Br V^^i) ^ (C7) 



G 



' eBVl^i 



f cos V — r{9 — UJrn) ( 1 + 



1 



1 + e cos V 



smt; 



ln{l + Br'V^^i] 



(C8) 



I have shown in paper I that at resonance, the numerical value of the logarithmic term 
In (1 + -B r^ ^rez) becomes so small that one can use the approximation ln(l+C) ~ C, ^ << 1 
to simplify equations (C5)-(C8). Using this approximation and also substituting for Um, 0, f 
and Vrei by (paper I) 
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Um ^ a-^/^ (l + -ecosl) , (C9) 

e ~ a-^/2 (l + 2ecos/) , (CIO) 

r ~ ea-^/^sin/ . (Cll) 

and 

Vrei ^ a-i/2e(l - ^cos^l) , (C12) 

respectively, TZl^TZq^T^i and TZg can be written as in equations (42)-(45). 
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